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Abstract. The possibility of observing phenomena peculiar to long-range 
interactions, and more specifically in the so-called Quasi-Stationary State (QSS) 
regime is investigated within the framework of two devices, namely the Free- 
Electron Laser (FEL) and the Collective Atomic Recoil Laser (CARL). The QSS 
dynamics has been mostly studied using the Hamiltonian Mean-Field (HMF) toy 
model, demonstrating in particular the presence of first versus second order out-of- 
equilibrium phase transitions from magnetized to unmagnetized regimes. Here, we 
give evidence of the strong connections between the HMF model and the dynamics 
of the two mentioned devices, and we discuss the perspectives to observe some 
specific QSS features experimentally. In particular, a dynamical analog of the 
phase transition is present in the FEL and in the CARL in its conservative regime. 
Regarding the dissipative CARL, a formal link is established with the HMF model. 
For both FEL and CARL, calculations are performed with reference to existing 
experimental devices, namely the FERMI@Elettra FEL under construction at 
Sincrotrone Trieste (Italy) and the CARL system at LENS in Florence (Italy). 
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1. Introduction 

Long-range interactions have now been shown to be central in a wide range of scientific 
contexts such as astrophysics [1], hydrodynamics [2] or nuclear physics [3]. However, 
the possibilities of investigating the long-range features via dedicated experiments are 
more restricted: Non- neutral plasmas [4 , cold atom and wave-particle systems jSj are 
among the most serious candidates. The purpose of this paper is to investigate the 
possibility of using existing set-ups based on the wave-particle interactions to probe 
long-range features of the dynamics, in particular out-of-equilibrium transitions. 

As an introduction to long-range interactions, let us start from the Hamiltonian 
Mean-Field (HMF) model J6j, a paradigmatic system on which many theoretical 
studies focused. This one-dimensional model describes the interaction of N particles 
on a circle through a collective field, which depends only on their phase 6j. This 
A^-body dynamics is described by the following Hamiltonian: 

N / 2 N \ 

j=i \ fe=i / 

associated to the canonical bracket in {9j,pj). Here, e = ±1 corresponds either 
to a ferromagnetic or an antiferromagnetic (-) system. In this model, the 
particles are collectively interacting through the so-called magnetization M = A/e"^ = 
(Sj ) /-^' since the dynamics of a single particle is given by: 

0j + eM siniOj - (jj) ^ 0. (2) 

Long-range systems can exhibit interesting equilibrium features, such as ensemble 
inequivalence (see e.g. [7] for the antiferromagnetic two-dimensional version of the 
HMF model or [8] for a recent review). However, the HMF model mainly revealed 
itself as a perfect playground to study out-of-equilibrium long-range features. Indeed, 
starting from generic non-stationary initial conditions, the system will typically have 
a fast transient dynamics until a nearly-stationary state, generally called Quasi- 
Stationary State (QSS), is reached: Not only this QSS dynamics substantially differs 
from the equilibrium one, but the system actually stays trapped in it for very long 
times [6j. 

More specifically, several authors actually demonstrated that the lifetime of 
the said QSS diverges when the number of particles in interaction increases. For 
example, numerical works report that the time of relaxation to equilibrium for the 
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Hamiltonian Mean- Field model scales as iV^-^ [5], in a regime of parameters yielding 
homogeneous QSS. To gain insight into the emergence of QSS, one can resort to a 
continuous picture, formal limit of the governing discrete Hamiltonian. A rigorous 
mathematical procedure leads to the Vlasov equation for the evolution of the single 
particle distribution function, the continuous representation of the particles density 
in phase space which is recovered when making the number of bodies N tend to 
infinity. The stability of QSS in the infinite N limit suggests that these latter 
states can be potentially interpreted as Vlasov stationary states, an ansatz that 
opens up the perspective for further analytical progress, a fact on which we shall 
return in the following. Operating in this context and explicitly accounting for 
finite size corrections beyond the idealized Vlasov picture, the authors of [10] proved 
rigorously that the relaxation of the A^-body system towards its deputed equilibrium, 
as driven by microscopic collision effects, would occur on time scales larger than TV, 
in qualitative agreement with the numerical evidences commented above. Clearly, 
QSS are supposedly the only regimes which are made experimentally accessible, in all 
physical situations where a large number of microscopic constituents evolve in mutual 
interaction. The experimental time of observation is in fact limited, and not sufficient 
to allow for equilibration. In this perspective, to unravel the puzzle of QSS and so 
build a comprehensive dynamical picture for their existence and evolution, represents 
a major challenge, with undoubtedly many practical implications. 

An important step forward explaining the presence of the QSS was eventually 
attained thanks to the theory of violent relaxation of Lynden-Bell (LB) [11]. This is 
a statistical theory which embeds self-consistently knowledge of the governing Vlasov 
dynamics. The approach is justified from first principles and allows to resolve the 
intermediate regime of the discrete A^-body evolution, when the system is presumably 
assimilable to a continuum Vlasov model, before finite size corrections come eventually 
into play. The theory is based on the maximization of the following entropic functional 
of the distribution function (DF) /: 



where /o describes the initial state of the system, whereas / stands for a coarse-grained 
distribution function of the final state, that one wishes to recover via a predictive 
approach. The above formulation holds for a two-step initial distribution function 
(water bag): / at time is equal to either zero or /g. Whereas the exact evolution 
according to the Vlasov equation imposes that the DF is only allowed to take and 
/o values at all times, the coarse-grained point-of-view implies a continuous DF / that 
is expected to be valid if one averages over small patches of phase space. As a side 
comment we notice that the functional ^ can be readily generalized to account for a 
continuous collection of different density levels, beyond the water-bag hypothesis. 

The maximization of s is performed under the macroscopic constraints of 
normalization, energy and momentum which are conserved by the dynamics. An 
underlying hypothesis to the theory is that the system explores in an ergodic-like 
fashion all states allowed by the constraints. The dynamical evolution of the Vlasov 
equation departs from that of a system sampling the equilibrium microcanonical 
ensemble, giving rise to different predictions which reflect the out-of-equilibrium nature 
of the problem. The application of the above predictive strategy to the study of 
the QSS dynamics of respectively the HMF model [12], free-electron lasers [13] and 
gravitational systems [2] has confirmed its adequacy. 




(3) 
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The LB approach also brought sorae new insights into the HMF phenomenology. 
For example, the abrupt change in the QSS magnetization when smoothly tuning the 
initial state of the system was interpreted as an out-of- equilibrium phase transition, 
which not only depends on the energy of the system - as it happens at equilibrium 
- but also on the precise way the system is prepared. More specifically, the initial 
magnetization of the system was shown to have an important role, and phase 
transitions of both first and second order could be observed depending on the value 
of this latter parameter. This memory effect - the system keeps track of the detail of 
its initial state for very long times - makes the QSS dynamics significantly richer than 
the equilibrium one. For example, as regards the HMF model, an out-of-equilibrium 
tricritical point was identified, which does not exist at equilibrium. 

The purpose of this paper is to determine whether some of the QSS features 
predicted for the HMF toy-model can be observed in experiments run for a dedicated 
class of devices. Motivated by this working hypothesis, we shall turn to considering 
the wide field of wave-particle interaction and focus in particular onto two different 
experiments, namely the Free-Electron Laser (FEL) and the Collective Atomic Recoil 
Laser (CARL). In both cases the dynamics reflects the long-range nature of the 
interaction, along the lines depicted above with reference to the simplified HMF 
setting. Operating in this framework, we will show that some features of the QSS 
dynamics, as those previously outlined, may be observed in direct experiments. 
Moreover, such properties though peculiar to the considered wave-particle dynamics, 
bear some reminiscent traits of the HMF model, to which both FEL and CARL 
are intimately connected. Eventually, the associated experimental set-ups are briefiy 
detailed, based on existing machines and current technology. 

Section [2] is devoted to FELs. The aim of such devices is to produce high-power 
short-wavelength light pulses by exploiting the radiation emitted by ultra-relativistic 
electrons when passing through the static and periodic magnetic field generated by 
an undulator. Starting from generic initial conditions, the wave power grows to a 
maximum, and then starts oscillating, keeping a lively exchange of energy with the 
particles, over times diverging with the number of particles, a characteristics of the 
QSS. As for the case of HMF, the QSS of a FEL depends not only on the energy 
of the system, but also on the details of its initial state. Thus, after presenting the 
FERMI@Elettra FEL, we discuss how to manipulate the electron beam to produce the 
sought different initial states. Finally, the dynamical transition present in the system 
is described. 

Section [3] is dedicated to discussing the Collective Atomic Recoil Laser (CARL), 
an experiment where a probe wave is amplified thanks to a grating of cold atoms 
(back)scattering photons of an incident pump laser beam. As for the FEL, its dynamics 
is dominated by long-range effects in the one-dimensional limit, an approximation 
which holds for the CARL experiment based at the European Laboratory for Non- 
linear Spectroscopy (LENS). We then focus on the conservative regime, when the 
dynamics formally reduces to that of the FEL: The possibilities to observe for the 
CARL the QSS phenomenology as depicted for the FEL is investigated. On the 
other hand, when the wave amplification takes place in a cavity, a damping has to 
be accomodated for: A formal link between this operational regime of CARL and 
the HMF dynamics is drawn, as well as the experimental perspectives to detect the 
associated out-of-equilibrium transitions. 

Finally, in Section |4l we discuss the measurements that could be performed 
for both CARL and FEL in order to unravel the imprint of QSS that indirectly 
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materializes in the existence of distinct out-of-equilibrium regimes. 

2. The Free-Electron Laser as a long-range interacting system 

FELs are powerful light sources able to deliver coherent pulses of photons over a large 
and tunable wavelength range. To that aim, ultra-relativistic electrons of energy 7 are 
injected into the periodic magnetic field (of period Xw and deflection parameter K) 
produced by an undulator, where they start to wiggle and emit synchrotron radiation 
around the following wavelength: 

+ (4) 

The light produced by the electrons traps the electrons themselves, resulting in 
a periodic modulation of the electrons' density (see FigH]) called bunching: This 
bunching is the source of the coherent emission. Eventually, under a resonant condition 
between the electrons and the wave, the strong interplay between coherent emission 
and particle trapping inside the wave potential leads to the nonlinear growth of the 
wave (see Fig[T]) and to the emission of a powerful light pulse. Due to the high energy 



1.5 
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Figure 1. Left: Electron phase-space in the QSS regime (at t = 50). Right: 
Normalized laser intensity versus normalized time t. Simulations performed with 
A'^ = 8000 particles, starting with a waterbag with bo = 0, Ap = 0.05, and a 
negligible intensity /q = 10~®. 



of the electrons (of order 1 GeV typically), the system can be in first approximation 
considered as one-dimensional, since the angle of the cone of light radiated goes as the 
inverse of the electrons energy. As for the radiation, it can generally be described by 
a mean-field wave, leading to the set of equations |15) : 




— = i V e^'^J + iSA 

where 9j is the phase of electron j with respect to the ponderomotive potential, pj its 
normalized energy, whereas A stands for the complex amplitude of the synchrotron 
radiation. The normalized variables are defined as 9j = {k + k^)zj — ujt ~ 5t, 
with Zj the position of particle j along the propagation axis, pj — (7^ — 7o)/p7oj 
7o the average electron energy, k and uj the radiation wavenumber and frequency, 
S = (70 ~ 1r)/p1q the detuning parameter and 7/j the resonant energy defined 
by Eq.(|4]). p — (///A)^/^(Au,au,/27rCT)^/^/27o is the so-called Pierce parameter. 
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= \/2K , I — rieiTTa'^ec the electron current, the electron density and Ia ~ 17kA 
the Alfven current. A corresponds to the normalized electric field of the wave, 
according to A = eo/{mc^'^oneeap), while the rescaled time t is given hyt — 2fc^pz, 
with z the position along the propagation axis. 

Following the HMF approach, we focus on the waterbag initial conditions, since 
they are a good description of the electron bunch as a first approximation [13]: the 
initial wave is initially of zero amplitude, while the particles are bunched into an 
homogeneous rectangle in the {9,p) phase-space, i.e. spread between —A9 and +A9 
in phases, and between — Ap and +Ap in momenta. Experimentally, shaping the initial 
electron bunch is part of the High Gain Harmonic Generation (HGHG) scheme [34] , 
where use is made of two distinct stages of interaction (see Fig(2|): In the first undulator 
sections, called the "modulator", electrons interact with an external coherent light 
source, e.g., a high-power laser (called the "seed"). Such an interaction induces the 
electron bunching at the seed wavelength, A^eed, and at the harmonics of the latter. 
In a second undulator section, tuned at one of the seed harmonics wavelength and 
called the "radiator" , electrons emit coherently. 



dispersive 

modulator X section radiator 




Figure 2. Schematic layout of the HGHG scheme: The electron beam is 
synchronized with the seed laser, which creates an energy modulation in the 
former inside the modulator. In the dispersive section, the energy modulation 
is converted into a spatial one, called micro-bunching. In the radiator, the micro- 
bunched electron beam emits coherently. 

In general, the wave grows, first quadratically, then exponentially, until it reaches 
a maximum and starts oscillating around an average value /. As for the electrons, 
they bunch together, thus allowing the coherent emission; the transfer of energy to 
the wave leads to a decrease in the electrons energy, spoiling the resonant condition 
([^. The amplification process stops when the particles are not any more in resonance 
with the wave (saturation) . Note however that when the energy spread Ap is too large 
(typically Ap > 1.5), the interplay between the wave and the electrons will not even 
trigger, and the wave amplification will not happen. 

The FERMIOELETTRA is a new EEL, presently under construction at the 
Sincrotrone Trieste laboratory. It aims at producing GW optical pulses in the 
10 — lOOnm range, thanks to the HGHG process: In this scheme, when considering 
a jump from Xseed in the modulator to its nth harmonic {Xrad — Xrnod/n, with n an 
integer) in the radiator, the bunching at the entrance of the radiator, is given by ^3] 

|fe„| = I < e"'' > I = e-^"'-''^V„(ndA7), (6) 

where <> is the average over the particles, J„ the n-th Bessel function of the first 
kind, cr-y the initial energy spread of the electron beam, and A7 the "coherent" energy 
spread generated by the modulation (see [33] for details, and Tab|T]for the FERMI 



CONTENTS 



7 



parameters). Here, d is the strength of the dispersive section, whose role is to convert 
the energy modulation into a spatial one. 

From now on, let us consider a HGHG configuration where the seed wave is 
Amod = 200nTO, associated to an harmonic jump of n = 2 {Xrad = lOOnm), with 
electron injected at 70 = 1760 into a. z = 18.4m-long radiator (see TabU] for the 
FERMI FEL parameters). Then, the maximum bunching is reached when the Bessel 
function J2 is maximized, that is for ndAj 3.05, and setting the n^a'^d'^ term close 
to zero. Thus, a dispersive section strength of d = 0.63/(ncr-y) leads to a decrease 
of 20% in \bn\ for the exponential term, and corresponds to a A7 w 3.5ct^. The 

initial and coherent energy spreads accumulating as a^^tot = \J <^'i^ + (^7)^/2, we get 

'^■y,tot ~ 0.12, whereas the bunching factor created is « 0.4. 



Table 1. Main parameters of FERMI modulator and radiator sections. 



Section 


Lw 


K 


70 


P 


A 




Modulator 
Radiator 


3m 
13.8-18.4m 


1-5 
1-5 


1760-2940 
1760-2940 


- 3.10-3 

- 3.10-3 


800-lOOnm 
100-20 nm 


0.035 
0.035-0.5 



As for the detuning, it is induced by shifting the resonant energy in the radiator 
from the average electron energy 70, according to the relation: 

6=1^^^. (7) 
P7o 

The / > saturated regime of the FEL was shown to be accurately described by 
the LB approach |T3]: The intensity and bunching reached by the laser are in good 
agreement with those determined by the maximization of entropy principle. It seems 
however not to apply to the non-lasing regime, when the resonance between the wave 
and the electrons is not satisfied anymore. A recent work [32] reported that regarding 
the LB principle, two solutions of the maximization problem exist: the one associated 
to a positive laser intensity is always entropically favored, but the system dynamics 
can actually be trapped in the vicinity of a zero-intensity solution, where the electrons 
stay unbunched. 

This dynamical trapping exhibits striking similarities with the phase transition 
encountered in HMF. Indeed, when monitoring control parameters such as the initial 
bunching 60 or the energy spread Ap, a transition from a J > regime to a / « 
one occurs, which can be sharp or smooth (see [32]). Let us focus on the parameters 
60 and detuning S: Fig|3] shows how the / > regime may abruptly end for bo = 0.1 
or values of the detuning 5 = 2 (top panels), whereas the transition to low-/ regimes 
is smooth for Bq = 0.5 01 S = 2.6 (bottom panels). It is worth noting that working 
at finite undulator length does not modify substantially the transition characteristics 
(dashed lines of Fig 121). 

The regions of parameters where each transition will occur is summarized in FiglJ] 
(left panel) , where the diagram of saturated intensity is depicted as a function of both 
parameters: It reveals that for bo below some be ~ 0.3, as well as for S below 6c ~ 2.3 
(top panels), the transition is sharp, whereas beyond these values, it turns out to be 
smooth. Thus, {be, Sc) represents the critical values of parameters beyond which the 
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b =0.1 CARL [LENSJ 

n 'i''. ' FEL[FERM1J 

QSS regime 





Figure 3. Intensity radiated versus control parameter {5 on the left [at fixed bg], 
and bo on the right [at fixed 5]). The intensity for the FEL (dashed lines) are 
simulated using the FERMI parameters (70 = 1760, A = lOOnm, p = 4.35.10"^, 
z = 18.4m, cTj^tot = 0.12); for the CARL (plain lines), the LENS parameters have 
been used {pc = 1000, t = l/is and energy spread ap = 20hk). The saturated 
intensity / of the QSS regime (dotted lines) is calculated as the average of the 
intensity between t = 50 and 100. Simulations performed with A'^ = 10000, and 
no initial wave. In the last picture, the CARL curve is represented three times 
larger. 



transition from lasing to non-lasing turns into a smooth one: Although it is in no way 
supported by an entropic approach as it was for HMF, it can be seen as a dynamical 
version of the tricritical point present in the latter toy-model. 

Regarding the possibilities to observe this peculiar phenomenology on the 
FERMI@Elettra FEL, a similar phase diagram has been plotted on Fig|3] (middle 
panel), which accounts for the finite interaction length of the machine. It reveals that, 
despite the deep saturation may not be well established within the undulator length 
available, the two areas where a sharp vs. smooth transitions could be observed are 
quite well separated. 

3. The Collective Atomic Recoil Laser 

CARL consists of a collection of cold two-level atoms driven by a far-detuned laser 
pump of frequency ujp which radiates at the frequency w ~ in the direction opposite 
to the pump [16j. In both the FEL and CARL systems the radiation process arises 
from a collective instability which originates from a symmetry breaking in the spatial 
distribution, i.e. a self-bunching of particles which group in regions smaller than the 
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Figure 4. Intensity as a function of {bo,S) in the QSS regime (left), and at the 
exit of the FERMIQElettra FEL (center) and CARL of LENS (right). Bright 
areas correspond to high intensities, dark ones to low-intensity regimes. The 
white lines represent the transition between the / > regime and the / one, 
with a transition either sharp (plain line) or smooth (dashed line). The white dot 
stands for the tricritical point {bc,Sc), when the transition in the intensity goes 
from sharp to smooth. Simulations performed in the same conditions as Fig[3] 



wavelength. In the limit in which the radiation pressure due to the pump laser can 
be neglected (for instance by largely detuning the pump frequency from the atomic 
resonance), the CARL is described by the same dimensionless FEL equations ([5|) 

Pj (8) 



(Ae*^^- + A*e-'^^) (9) 



^3 

di 
dt 

f =l;^e--^+.M-.A (10) 
1=1 

however with the presence of a damping term — Kj4 in the field equation, accounting 
for radiation losses from a ring cavity surrounding the atoms. Although CARL 
and FELs evolve with a similar dynamics, the dimensionless variables of the two 
systems, and consequently also the typical timescales, are very different. In CARL, 
the phase and the normalized momentum of the atoms j are 9j — 2k{zj(t) — {vz)ot) 
and Pj = m{vzj{t) — {vz)o)/{2hkpc) (where Zj{t) and Vzjit) are the position and 
velocity of the jth atom along the direction of the scattered field and {vz)o is the 
average initial velocity), whereas A stands for the normalized complex amplitude 
of the radiation field, A = {eQ/hujUaPcY^'^EQ, where ria is the atomic density. 
The scaled time \s i = {SuJrecPc)^, where uJrec = hk^ /2ra is the recoil frequency, 
5 = (lu — LOp — 2k{vz)o) / (SujrecPc) IS the pump-probe detuning, k = Kc/ (SuJrecPc) is 
the scaled loss of a ring cavity with length Lcav, transmission T and Kc = cT/Lcav, 
and finally pc — (T /8){c(7Qna/ A'^ojf^^Y^'^ {I / IsatY^^ , where F is the natural decay rate 
if the excited state, (Tq = 3A^/27r is the scattering cross section, A = ujq — ujp is the 
pump-atom detuning, / is the pump intensity and Isat — ?iwF/2ctq is the saturation 
intensity. 

There have been different experiments that have observed the CARL effect in 
room temperature gases [T7], cold atomic samples from Magneto-Optical traps (MOT) 
[H m] or Bose- Einstein condensates (BEC) [IHl IHl US- CARL experiments were 
performed either in high-finesse optical cavities ITHl or in free space [H], where the 
effect was originally interpreted as Superradiant Rayleigh Scattering [20]. However, 
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matter wave grating 




Figure 5. Schematics of superradiant light scattering from a Bose-Einstein 
condensate. An elongated BEC is illuminated by a far off-resonant laser beam 
(pump beam) with frequency uj and wavevector k directed along its axial direction. 
After backscattering of photons with kg ~ —k and the subsequent recoil of atoms, 
a matter wave grating forms, due to the quantum interference between the two 
momentum components of the wavefunction of the condensate. The effect of this 
grating is to further scatter the incident light in a self-amplifying process. 

it has later been empliasized that these experiments in free space can be seen as a 
CARL process in the superradiant regime [321 HI] ■ 

At LENS, the CARL experiment [IS] is reahzed with a cigar-shaped BEC of ^''Rb 
produced in a lofFe-Pritchard magnetic trap by means of RF-induced evaporative 
cooling. After 2 ms of free expansion, when the magnetic trap field is completely 
switched off and the atomic cloud still has an elongated shape (at this time the radial 
and axial sizes of the condensate are typically 10 and 70 /xm, respectively), a square 
pulse of light is applied along the z-axis (see FiglS]). The size of the laser beams is 
larger than 0.5 mm, far larger than the condensate free fall during the interaction with 
light. In this geometry the CARL process causes the pump light to be backscattered 
and the self- amplified matter- wave propagates in the same direction as the incident 
light. 

This experiment allows for a great flexibility in the preparation of the initial state 
of the system: for example, it is possible to prepare the atoms in an initially bunched 
state by imposing an electromagnetic standing wave before the pump laser is activated. 
Regarding the momenta spread, it can be varied by cooling only partially the atoms. 
Finally, a non-zero pump-probe detuning d can be induced by giving the atoms an 
initial momentum, which can be up to 1000fi.fc: Indeed to change the initial detuning 
of the probe, one can either change the initial atomic momentum or change the cavity 
frequency since conservation of energy imposes the frequency of the backscattered 
photon [26] . 



Table 2. Main parameters of LENS CARL. 



Wrec PC 


t 


S K Ap 


5000Hz 100-1000 


l/is-5ms 


-5 - 5 0.025-10^ 0.01-5 



Equilibrium dynamical phase transitions were studied in the CARL context under 
different hypotheses: For example, using a thermalization hypothesis for the atoms, 
the CARL will reach equilibria where (equilibria) phase transitions can occur |27j . 
In the case where dissipation is counterbalanced by a stochastic Langevin force |19) . 
steady states can be identified, as well as a phase transition analogous to the one that 
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occurs in the Kuramoto model [35], yet with a self-generated collective oscillation 
frequency for the CARL. Both models were also connected when the strong friction 
is balanced by diffusion a regime for which steady states and phase transitions 
were also predicted [5U] . Yet, the equilibrium condition remains a condition which 
is hard to satisfy experimentally. In this context, the study of the Quasi- Stationary 
regime, which are not iV-body equilibrium regimes but rather stationary states of the 
associated Vlasov dynamics, is all the more relevant experimentally. 

Let's first consider the configuration where the CARL amplification is realized 
inside a ring cavity: the pump and probe light fields are counterpropagating modes 
of the ring cavity and the interaction time of the light fields with the atoms can 
be enhanced by several orders of magnitude, which supports the amplification. 
Consequently, most of the CARL experiments carried out up to date employed ring 
cavities |18[ I19[ [ ^ . Furthermore, in new experiments, especially in the microtrap- 
based ones, it appeared to be very profitable to introduce a high finesse optical cavity 
inside the vacuum system: Indeed, in this environment, Fabry-Perot cavities have 
already been demonstrated to yield very high finesses |31j . 

Back to the dynamical equations (f8l-fT0|). it is clear that if the dissipation term k is 
small {k ^ 1), the CARL is governed by the same equations as the FEL. In particular, 
when monitoring the initial bunching and the detuning parameter, it should exhibit 
a dynamical transition similar to the one described in SecH] Simulations realized for 
a good cavity (k = 0.025) reveal that the sharp/smooth transitions could indeed be 
observed in the relevant range of parameters (see FigH]). 

On the other hand, the CARL experiment can be realized either without cavity, 
or within a low-quality one: This is at the cost of a high damping of the wave, 
modelled by the k > term in Eq. pTI| . In this case, if the dissipation is large enough 
(k '-^ 1), an adiabatic treatment of the wave dynamics can be performed [23[|25| . which 
corresponds to setting dA/dt = 0, which yields 

A^^—-ye-''^. (11) 

K-iSN^ ^ ' 

j 

Plugging this expression into the particles equations of motion ^ , we get the following 
equation for the momentum of particle j (the equation for its position is unchanged) : 

2 1 ^ 

Pi = s'^ + k'^ N^^^ -e,n)- K cos {9 j~e,n)). (12) 

m 

If we now consider the extra limit where the detuning S is large with respect to k 
(I (5 1 ^ k), Eq. ([T2|) simply turns into 

^^^■ = ^^E«in(0,-0„o. (13) 

Then, using the following normalization 



^AS\/2, 



(14) 



i = ty/\S\/2, H = \S\H/2, 

the system can be mapped into the HMF model ([T]). Anew, e = —sign{S) positive 
(negative) describes an (anti)ferromagnetic interaction. Numerical simulations 
confirm that in the above mentioned limit, despite small differences in the dynamics, 
the CARL and HMF yield similar QSS regimes (see FigH]). 

This formally links the CARL dynamics to the HMF model, in either its 
ferromagnetic or antiferromagnetic form, making in particular the bridge between 
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Figure 6. Comparison between the dynamics of the CARL (plain lines) and 
HMF (dashed lines) dynamics, in magnetized regime (left: 60 = 0, Ap = 0.1) and 
unmagnetizcd regime (right: bo = 0-94, Ap = 1.5). Simulations performed with 
N = 10000 particles; CARL parameters: k = 1 and S = —4. 



the rich phenomenology predicted for the HMF model and a possible experimental 
realization. For example, regarding the ferromagnetic case, FiglT] depicts the phase 
diagram of HMF in the (60 , U) plane as predicted by the LB prescription, and those 
obtained by direct A^-body simulations of both the CARL dynamics and the HMF 
model. They reveal an excellent agreement between the CARL phase diagram and 
its reduced counterpart, the HMF, as well as with the Lynden-Bell approach. Note 
that the presence of fringes in the large U part in the diagrams are due to the short 
time considered: yet, although the QSS regime is not well established, the transition 
line is already present and in good agreement with the LB prescription. This allows 
to conclude on the presence of an out-of-equilibrium phase transition in the CARL 
device, such as predicted for the HMF model [35] . 



0. 



0. 




0.2 0.4 0.6 0.8 0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



^^0 ^0 

Figure 7. Bunching factor as a function of the initial bunching and energy of 
the system, as predicted by Lynden-Bell's violent relaxation approach (left), and 
by Af-body simulations of the CARL (center) and HMF (right) dynamics at finite 
length. The white line stands for the transition as predicted by Lynden-Bell, of 
the first order type below Mc ~ 0.17, and of the second order kind above. Af-body 
simulations realized with N = 10000 until t = 40; CARL simulations performed 
with K = 0.5 and 5 = —5. Note that the energy U here refers to the normalized 
energy U = Ap^/e + (1 - b^)/2. 



4. On the possibility to observe the QSS signature 



In this paper, we investigated the out-of-equilibrium dynamics of two long-range 
models which admit corresponding experimental implementations. We shall now 
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conclude the discussion by elaborating on the experimental possibilities to verify 
the correctness of the proposed picture and in particular detect the predicted phase 
transitions. 

As concerns the FEL, the particles phase-space is not experimentally accessible 
since the electrons have relativistic energies. However, by monitoring the laser 
intensity it should be in principle possible to detect the transition, as demonstrated 
in SeclSl Depending on the specific initial condition, ultimately characterized by an 
assigned electron bunching amount, the laser can gain in potency or, alternatively, 
keep its off mode. Experimentwise, by setting the initial bunching at either zero (no 
interaction in the modulator) or high values (60 ~ 0-4), and tuning the electrons 
energy, one should in principle observe either a sharp or a smooth transition, as 
indicated by the simulations results. This is an interesting feature, indirect signature 
of the QSS existence, which bears an intringuing similarity with the HMF behaviour. 
In this latter case, however, the presence of an out-of-equilibrium phase transition of 
both first and second order types was demonstrated on solid theoretical grounds. A 
similar theoretical justification is still lacking with reference to the FEL model. 

The CARL device, in the regime of small dissipation, is predicted to exhibit 
a transition, similar to that displayed by the FEL. Again, and in analogy with the 
above, we suggest that the presence of the transition could be successfully evidenced by 
recording the radiated intensity, under different experimental conditions. In addition, 
the CARL could also allow to access direct information on the particles dynamics in 
phase space: First, the atoms density can be recorded, and the presence of fringes 
used to quantify the degree of bunching. Then, letting the atoms expand after the 
CARL process ends, and recording their late positions, one can resolve the momentum 
distribution: Such a diagnostic would translate into an independent tool to bring into 
evidence the transition between distinct dynamical regimes. 

Finally, it is when operating the CARL device in the dissipative HMF-like regime 
that the possibility to measure the phase-space structures could be of paramount 
importance. The fringes, which means bunching, are in principle strongly correlated 
to the radiated power (see Ea lTT|) . By accessing the distribution of momenta one could 
eventually detect the two bumps that are predicted to occur in presence of unbunched 
QSS, and which are believed to correspond to counter-propagating clusters of particles. 
More generally, and with reference to the ferromagnetic case [121 ISS] we shall be 
interested in accurately investigating the overall velocity profile whose characteristics 
have been object of vigorous debates [37l |38] . 
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